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Abstract — A  new  integrated  collector  storage  (ICS)  concept  for  low-temperature  solar  heating  of  water  is 
described.  The  solar  energy  is  stored  in  a  salt-hydrate  phase-change  material  (PCM)  held  in  the  collector 
and  is  discharged  to  cold  water  flowing  through  a  surface  heat  exchanger  located  in  a  layer  of  stationary 
heat  transfer  liquid  (SHTL),  floating  over  an  immiscible  layer  of  PCM.  A  theoretical  model  for  the 
charging  process  of  the  proposed  integrated  collector  is  presented.  The  model  assumes  one-dimensional 
transient  heat  conduction  in  the  PCM  and  SHTL  layers  and  neglects  the  effect  of  convection  heat  transfer 
in  these  regions.  The  model  was  solved  numerically  by  an  enthalpy-based  finite  differences  method  and 
validated  against  experimental  data.  The  results  of  parametric  studies  on  the  effect  of  the  transition 
temperature  and  of  the  thickness  layer  of  the  salt-hydrate  PCM  on  the  thermal  performance  of  the  charging 
process  are  also  presented. 


1.  INTRODUCTION 

The  integrated  collector  storage  (ICS)  concept 
is  considered  to  be  a  promising  direction  for 
increasing  the  economic  feasibility  of  low- 
temperature  solar  systems  for  heating  water  for 
domestic,  agricultural  and  industrial  applica¬ 
tions.  A  system  of  this  type  combines  collection 
and  storage  of  thermal  energy  in  a  single 
unit.  Compared  with  the  conventional  domestic 
water  heating  system,  the  integrated  collector 
has  the  advantage  of  simplicity,  both  in  erection 
and  in  operation.  However,  it  has  a  relatively 
low  efficiency  (Bar-Cohen,  1978;  Shmidt  et  al., 
1988).  Boy  et  al.  ( 1987)  proposed  an  ICS  system 
based  on  a  salt-hydrate  phase-change  material 
(PCM)  as  an  appliance  for  providing  hot  water 
instantaneously.  They  demonstrated  that  the 
thermal  efficiency  of  such  system  could  be 
improved  significantly  by  incorporating  an 
appropriate  PCM  device.  However,  in  their 
system  the  salt-hydrate  PCM  was  encapsulated 
in  a  special  corrugated-fin  heat  exchanger,  which 
increased  the  cost  of  the  system. 

In  the  present  work  we  propose  a  new  type 
of  integrated  collector  module  that  does  not 
require  a  water/PCM  surface  heat  exchanger 
but  uses  a  water/organic  liquid  (oil)  heat 
exchanger  instead.  A  schematic  representation 
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of  the  proposed  apparatus — the  ICS-PCM 
module — is  given  in  Fig.  1.  The  system  consists 
of  a  thermally  insulated  plastic  container,  con¬ 
taining  a  layer  of  stationary  heat  transfer  liquid 
(SHTL),  floating  over  a  layer  of  an  immiscible 
salt-hydrate  PCM  (the  density  of  the  SHTL  is 
obviously  lower  than  that  of  the  liquid  salt 
hydrate).  A  finned  heat  exchanger  located  in  the 
SHTL  is  used  to  heat  cold  water  during  the 
discharging  process.  The  SHTL  is  covered  with 
a  metal  plate  coated  with  a  selective  absorbing 
film  and  with  a  transparent  insulating  lid.  As 
the  PCM  melts,  its  density  decreases  and  its 
volume  thus  increases.  The  SHTL  displaced  as 
a  consequence  of  this  expansion  is  stored  in  a 
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small  interconnecting  chamber,  designated  the 
expansion  vessel.  A  removable  cover  made  of  a 
material  with  a  low  thermal  conductivity  and 
coated  on  one  side  with  a  reflective  material 
may  be  used  as  an  optional  extra  to  increase 
solar  energy  absorption  during  the  day,  as 
shown  in  Fig.  1.  Tt  may  also  be  used  as  addi¬ 
tional  cover  to  reduce  heat  loss  during  the  night. 

The  charging  process  taking  place  during  the 
day  is  based  on  the  absorption  of  the  solar 
radiation  by  the  selective  surface  and  heat 
transfer  to  the  PCM  through  the  SHTL  layer, 
mainly  by  conduction.  This  process  causes  the 
PCM  to  melt.  The  discharging  process,  which 
is  based  on  the  falling  crystal  method  (FCM) 
(Greene,  1978;  MacCraken,  1981;  Korin  and 
Mikic,  1989)  takes  place  as  follows.  As  the 
temperature  at  the  SHTL/PCM  interface  drops 
below  the  solidification  temperature,  solid 
species  are  formed  in  the  PCM.  Under  appro¬ 
priate  design  and  operation  conditions,  the  solid 
crystal  seperate  from  the  SHTL/PCM  interface 
and  sink  to  the  bottom  of  the  container  under 
the  influence  of  gravity.  “Fresh”  liquid  species 
thus  come  into  contact  at  the  interface,  and  the 
cycle  is  repeated.  The  direct  contact  between 
the  SHTL  and  liquid  salt  hydrate  facilitates 
good  thermal  contact.  Heat  released  in  this  way 
is  transferred  to  the  cold  water  flowing  through 
a  finned-coil  heat  exchanger  located  in  the 
SHTL  layer. 

The  main  advantages  of  the  proposed  ICS- 
PCM  method  are:  (a)  the  system  does  not 
require  an  expensive  phase-change  surface  heat 
exchanger;  (b)  the  water/oil  heat  exchanger  does 
not  come  into  contact  with  the  salt-hydrate  and 
therefore  corrosion  problems  are  prevented; 
(c)  in  cold  areas  anti-freeze  protection  is  not 
required;  and  (d)  the  system  is  simple,  and  there 
are  no  special  requirements  for  maintenance 
and  operation.  The  main  disadvantages  of  the 
system  are:  (a)  application  of  this  method  is 
limited  to  those  salt  hydrate  PCMs  that  have 
large  differences  between  their  solid  and  liquid 
densities;  and  (b)  the  discharge  process  depends 
on  the  kinetics  of  crystallization  of  the  chosen 
salt  hydrate  PCM.  [Relevant  mass  transfer 
problems  with  their  possible  solutions  related 
to  the  discharging  process  are  given  in  Korin 
and  Mikic  (1989).]  These  drawbacks  limit  the 
possible  applications  of  this  method  to  special 
cases. 

We  thus  intend  to  develop  the  system  for 
solar  heating  of  water  to  relatively  low  temper¬ 
atures  for  applications  such  as  space  and  root 


zone  heating  in  protected  cultivation.  For  exam¬ 
ple  in  a  closed  greenhouse  located  in  an  arid 
region  (such  as  in  the  area  of  Beer-Sheva,  Israel), 
the  air  temperature  inside  the  building  may  rise 
to  >35°C  on  a  typical  clear  winter’s  day  and 
drop  during  the  night  to  values  close  to  the 
ambient  minimum  night  temperature,  which  is 
usually  in  the  range  of  4-5°C  (Korin  et  al., 
1987).  On  a  clear  cold  night,  due  to  intensive 
radiation  to  the  sky,  the  minimum  night  temper¬ 
ature  may  even  fall  below  0°C.  To  obtain  the 
desirable  indoor  temperature,  which  is  specific 
to  each  type  of  crop,  (e.g.  melons  12-35°C, 
tomatoes  10  30°C,  roses  16-28°C)  the  green¬ 
house  has  to  be  equipped  with  cooling  and 
heating  systems.  In  many  cases,  a  low-temper¬ 
ature  energy  source  (18-25°C)  would  be  suffi¬ 
cient  for  keeping  the  night  temperature  above 
the  required  minimum.  For  example,  to  main¬ 
tain  the  air  temperature  above  12°C  when  the 
outside  temperature  is  0-5°C,  a  heat  source  of 
20°C  would  be  sufficient.  The  ICS-PCM  system 
could  be  used  to  supply  warm  water  at  a 
relatively  low  temperature  for  keeping  the  night 
temperature  in  the  greenhouse  above  the  mini¬ 
mum  required.  The  system  could  also  be  used 
in  combination  with  an  other  solar  system,  such 
as  passive  solar  plastic  water  sleeves  (Zabeltitz, 
1989),  or  a  selective  radiation  filter  solution  in 
forced  circulation  (Kopel  et  al.,  1991). 

This  paper  focuses  mainly  on  the  study  of 
solar  charging  process  in  the  ICS-PCM  system, 
and  includes  the  following  aspects:  ( 1 )  a  theoret¬ 
ical  model  for  thermal  analysis  of  the  charging 
process  in  the  ICS-PCM;  (2)  validation  testing 
of  the  theoretical  model  versus  data  obtained 
from  an  experimental  module;  and  (3)  theore¬ 
tical  parametric  studies  related  to  the  phase- 
change  transition  temperature  and  the  thickness 
of  the  PCM  layer. 

2.  THEORETICAL  MODEL  FOR  THE 
CHARGING  PROCESS 

A  schematic  description  of  the  ICS-PCM 
during  the  energy  charging  period  is  presented 
in  Fig.  2.  The  model  relates  to  a  multilayer 
system  consisting  of  the  following  elements:  a 
transparent  cover,  an  air  layer,  a  selective 
absorber  plate,  a  SHTL,  a  PCM  salt-hydrate 
layer,  and  a  thermal  insulation  layer.  The  model 
assumes  temperature-dependent  thermophysical 
properties  and  time-dependent  meteorological 
parameters,  i.e.  solar  intensity,  ambient  temper¬ 
ature,  and  wind  velocity. 
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Fig.  2.  Schematic  presentation  of  the  theoretical  energy 
charging  model. 


Since  the  heat  capacities  of  the  transparent 
cover  and  the  air  layers  are  relatively  small  with 
respect  to  the  oil  and  the  PCM,  and  since  the 
time  constants  for  the  thermal  response  of  these 
layers  are  relatively  short  with  respect  to  the 
solar  cycle  period  (Shmidt  et  al.,  1988),  only  the 
thermal  resistance  to  heat  flow  through  these 
layers  are  taken  into  account  in  the  ICS-PCM 
mathematical  model.  The  temperatures  of  the 
transparent  cover  and  the  air  layer  are  consid¬ 
ered  to  be  uniform  and  time  dependent. 

Convection  heat  transfer  in  the  SHTL  and 
the  PCM  layers  is  neglected,  and  transient,  one¬ 
dimensional  heat  conduction  is  assumed.  For 
the  PCM  layer,  it  is  further  assumed  that  phase 
transition  in  the  salt  hydrate  occurs  over  a  range 
of  temperatures.  Thus,  the  PCM  layer  is  divided 
into  three  sublayers:  a  liquid  sublayer,  a  solid 
sublayer,  and  the  sublayer  undergoing  phase 
transition.  It  is  assumed  that  all  the  thermophys¬ 
ical  properties  of  the  PCM  are  constant  and 
uniform  in  each  sublayer.  However,  the  specific 
heat  in  the  phase  transition  sublayer  is  an 
effective  thermophysical  property  including  the 
latent  heat  effect  (Rabin  and  Korin,  1993). 

On  the  basis  of  the  above-described  assump¬ 
tions,  the  governing  equations  in  the  SHTL  and 
in  each  of  the  PCM  sublayers  are  similar  and 
may  be  presented  as: 


d2Tt  1  dTj 
dx2  <Xi  8t 


(1) 


where  i  is  the  layer/sublayer  index,  which  can 
be:  f  £,  m  or  s  for  the  SHTL,  liquid  PCM, 
melting  PCM,  or  solid  PCM,  respectively. 

The  absorber  plate  has  a  relatively  low  ther¬ 
mal  resistance  to  heat  flow  (metal)  and  a  very 
low  heat  capacity  (a  thin  plate)  with  respect  to 
the  other  components.  Thus,  it  is  assumed  that 
the  absorber  plate  acts  as  a  time-dependent 


plane  heat  source.  Neglecting  the  thermal  mass 
of  the  absorber  plate,  the  upper  boundary  condi¬ 
tion  for  eqn  (1)  could  be  obtained  from  energy 
balance  on  the  absorption  plate,  which  is  given 
by  eqn  (2) 


8Tt  -  Tt\x. o 
~dx+  Rc  +  R 


+  AaqT  =  0 


(2) 


where  the  parameter  a  is  an  effective  absorption 
coefficient  of  the  integrated  system  (measured 
experimentally),  R ^  is  the  convective  thermal 
resistance  between  the  surface  cover  and  the 
environment,  which  is  given  by: 


Roo 


1 

h^A 


(3) 


where  hx  is  the  convection  heat  transfer  coeffi¬ 
cient  (in  W  m-2  K”1)  determined  from  eqn  (4), 
in  which  U  is  the  wind  velocity  (in  ms”1) 
(Macadams,  1954): 


hx  =  6.2  +  1.4(7.  (4) 


Rc  is  the  equivalent  thermal  resistance  of  both 
the  transparent  cover  and  the  stationary  air 
layer  between  the  cover  and  the  absorber  plate 
and  is  given  by: 


k.-± 


(5) 


where  kc  is  an  effective  thermal  conductivity 
coefficient. 

The  inner  boundary  conditions  at  the  inter¬ 
faces  between  the  STHL  layer  and  the  PCM 
sublayers  present  continuity  in  temperature  and 
heat  fluxes: 


=  T]\*-h 

(6) 

dT.\ 

-  c.;i.  „ 

(7) 

where  the  index  i  may  represent  the  layers/ 
sublayers  /,  £  or  m,  and  the  index  j  represents 
£,  m  or  s,  respectively. 

The  lower  boundary  condition  is  presented 
as  a  thermal  insulation: 


The  initial  temperature  distribution  is  consid¬ 
ered  to  be  a  function  of  one  dimension  and  is 
represented  by: 

T(x,  f  =  0)  =  g(x)  0  <x<L.  (9) 

The  phase  change  problem  defined  above 
using  an  effective  specific  heat  is  based  on  the 
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enthalpy  approach.  Thus,  the  solution  of  this 
problem  was  performed  by  a  finite  difference 
and  enthalpy-based  numerical  scheme  (Rabin 
and  Korin,  1993).  This  numerical  solution  was 
found  to  be  an  efficient  numerical  technique 
for  a  mulidimensional  problem  under  various 
boundary  conditions.  The  proposed  solution 
provides  the  temperature  distribution  in  the 
system.  However,  the  locations  of  the  PCM 
interfaces  have  to  be  interpolated  from  the  tem¬ 
perature  distribution.  Furthermore,  it  should  be 
noted  that  there  can  be  one  or  more  interfaces 
at  the  initial  time. 

3.  EXPERIMENTAL  SETUP 

An  experimental  model  of  the  ICS-PCM 
system  was  designed  and  build  for  an  investiga¬ 
tion  of  the  concept  and  of  the  mathematical 
model  in  the  stage  of  the  charging  process.  The 
apparatus  was  designed  to  simulate  a  center 
part  of  the  ICS-PCM  prototype,  where  the  heat 
flow  is  one-dimensional  and  where  there  are  no 
end  effects.  Thus,  the  experimental  model  was 
designed  to  have  the  same  height  as  the  proto¬ 
type  system.  One-dimensional  heat  flow  in  the 
experimental  model  was  achieved  by  insulating 
the  model  with  a  thick  layer  of  a  material  hav¬ 
ing  a  very  low  thermal  conductivity  coefficient 
(polyurethane). 

A  schematic  cross-section  of  the  components 
of  the  experimental  model  are  shown  in  Fig.  3. 
The  experimental  system  consisted  of  a  10  mm 
thick  Plexiglas  container,  0.3  x  0.3  m  in  cross- 
section  and  0.131  m  in  height.  The  sides  of 
the  container  were  thermally  insulated  with  a 


40  mm  thick  layer  of  polyurethane.  The  con¬ 
tainer  was  filled  with  13.4  kg  of  a  salt  hydrate 
PCM,  composed  of  a  eutectic  mixture  (in  wt%) 
of:  48.0  CaCl2, 4.5  KC1, 0.4  NaCl  and  47.1  H20, 
with  the  addition  of  1  wt%  BaCl2  •  2H20  as  a 
nucleating  agent.  This  PCM  has  a  phase  trans¬ 
ition  temperature  interval  of  27-29°C,  a  latent 
heat  of  164.5  kJ/kg,  and  solid  and  liquid  densi¬ 
ties  of  1610  and  1490  kg/m3,  respectively  (Korin, 
1982;  Wolf  et  al,  1984;  Korin  et  al.,  1985). 
Paztherm  22  oil  (1460  cm3),  used  as  the  SHTL, 
was  poured  onto  the  PCM  in  the  solid  state. 

Before  each  experiment  the  system  was  kept 
in  a  constant  controlled  temperature  compart¬ 
ment  for  about  three  days  so  that  the  PCM 
could  solidify  completely.  Temperatures  profiles 
were  measured  by  copper-constantan  thermo¬ 
couples  (type  T)  mounted  on  three  3  mm  dia 
Plexiglas  bars,  each  containing  five  sensors.  The 
container  was  covered  with  a  metal  plate  coated 
with  a  selective  absorber  material.  A  10  mm 
thick  transparent  polycarbonate  panel  was 
placed  over  the  cover  to  provide  thermal  insula¬ 
tion.  The  oil  that  was  pushed  out  of  the  system 
during  the  melting  was  allowed  to  flow  through 
3  mm  copper  pipes  and  rubber  piping  into  a 
500  ml  graduated  vessel.  By  measuring  the 
volume  of  oil  that  flowed  out  of  the  system  the 
locations  of  the  liquid/solid  interfaces  were 
determined,  as  a  function  of  time,  for  further 
validation  of  the  model.  A  summary  of  the 
relevant  thermophysical  properties  of  the  com¬ 
ponents  of  the  experimental  ICS-PCM  system 
is  given  in  Table  1 . 

Meteorological  data  (ambient  temperature, 
solar  intensity,  and  wind  velocity)  were  obtained 


Fig.  3.  Schematic  cross-section  of  the  experimental  system. 
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Table  1.  Typical  thermophysical  properties  of  the  materials 
of  the  ICS  system 


Thermal  conductivity 

PCM  in  a  liquid  phase 

0.5t 

Wm-’r1 

PCM  in  a  solid  phase 

0.6f 

Plexiglas 

0.1 

Polyurethane 

0.04 

Absorber  plate 

43 

Air 

0.028 

(50°C) 

Transparent  cover 

0.042 

Paztherm  22  (oil) 

0.145 

Volumetric  specific  heat 

PCM  in  a  liquid  phase 

4. 16f 

MJ  m-3  K_1 

PCM  in  a  solid  phase 

2.88f 

Absorber  plate 

3.69 

Air 

0.0013 

Paztherm  22  (oil) 

1.62 

Density 

PCM  in  a  liquid  phase 

1600t 

kgm-3 

PCM  in  a  solid  phase 

1800t 

Paztherm  22  (oil) 

860 

Volumetric  latent  heat 
MJm-3 

PCM 

255t 

tAverage  values  of  salt  hydrates. 


from  a  meteorological  station  located  near  the 
experimental  area.  The  solar  energy  absorbed 
by  the  system  was  determined  by  a  measurement 
of  the  direct  and  diffuse  radiation  reflected  from 
the  upper  surface  of  the  system.  Uncertainty  of 
temperature  measurements  in  the  experimental 
apparatus  was  estimated  to  be  in  the  range 
of  ±0.5°C. 

4.  RESULTS  AND  DISCUSSION 
4.1.  Experimental  validation  tests 

Several  energy  charging  experiments  were 
conducted  in  an  arid  zone  area — Beer-Sheva, 
Israel  (31°15'N,  45°34'E,  315  m  MSL)— during 
the  wintertime.  Results  of  these  experiments 
were  used  for  validation  of  the  mathematical 
model  as  well  as  for  an  experimental  investiga¬ 
tion  of  the  ICS  concept. 

Typical  validation  test  results  of  the  mathe¬ 
matical  model,  from  one  of  these  experiments 
conducted  on  a  typical  winters  day,  are  given 
below.  Figure  4  presents  a  comparison  of  the 
one-dimensional  temperature  distribution,  as 
calculated  by  the  numerical  solution  of  the 
mathematical  model,  and  the  measured  data.  It 
can  be  seen  that  there  is  a  very  good  agreement 
between  the  numerical  results  and  measured 
data  in  the  PCM  sublayers.  The  numerical 
solution  in  the  SHTL  layer  is  less  accurate, 
especially  in  the  afternoon,  probably  due  to  the 
fact  that  the  mathematical  model  neglects  heat 
transfer  by  convection  in  this  liquid  layer.  The 
relatively  low  accuracy  of  the  numerical  solution 


in  the  SHTL  has  a  minor  influence  on  the 
thermal  analysis,  since  the  thermal  mass  of 
the  SHTL  is  considerably  less  than  that  of  the 
PCM.  On  the  other  hand,  the  high  accuracy  of 
the  numerical  solution  in  the  PCM  layer  is  of 
major  importance,  since  it  helps  us  in  estimating 
the  fraction  of  the  liquid  phase  in  the  PCM  and 
thus  in  estimating  the  amount  of  energy  stored 
in  the  system. 

Predicted  values  for  liquid/solid  interfaces 
and  the  energy  stored  in  the  system  were  in 
good  agreement  with  experimental  data,  when 
considered  for  engineering  design  purposes.  The 
results  show  that  on  a  typical  winters  day,  the 
salt-hydrate  started  melting  at  about  09:30,  and 
the  overall  daily  solar  absorption  efficiency 
(defined  as  the  ratio  between  the  total  energy 
absorbed  by  the  PCM  and  the  total  daily  solar 
radiation  reaching  the  collector  area)  was 
about  70%. 

4.2.  Parametric  studies 

The  phase  transition  temperature  interval  and 
the  thickness  of  the  salt-hydrate  PCM  layer  are 
important  design  parameters  of  the  ICS-PCM 
system.  In  the  determination  of  the  optimal 
values  of  these  parameters  for  a  specific  applica¬ 
tion  both  the  charging  and  the  discharging 
processes  have  to  be  considered.  In  the  present 
analysis,  the  validated  theoretical  model  was 
used  to  study  the  effect  of  these  parameters  on 
the  ICS-PCM  system  during  the  charging  pro¬ 
cess.  In  the  light  of  the  possible  application  of 
night  space  heating  in  greenhouses  in  an  arid 
zone,  the  study  of  the  phase  transition  parameter 
was  limited  to  a  range  between  15-35°C.  The 
study  was  performed  with  a  series  of  PCMs 
that  are  based  on  salt-hydrate  mixtures  of 
CaBr2*6H20  and  CaCl2*6H20.  The  unique 
features  of  these  PCMs  are:  (1)  the  various 
mixtures  solidify  as  solid  solutions  and  melt 
completely  over  a  narrow  temperature  range  of 
about  2°C  to  give  a  homogeneous  liquid  solu¬ 
tion;  (2)  with  suitable  nucleating  agents  the 
solidification/melting  process  takes  place  with¬ 
out  going  through  an  incongruent  phase  trans¬ 
ition  point,  so  that  the  phase  change  process  of 
these  mixtures  is  totally  reversible  (Korin,  1982); 
(3)  by  selecting  a  suitable  CaCl2-6H20  to 
CaBr2-  6H20  weight  ratio,  it  is  possible  to 
adjust  the  transition  temperature  within  a  range 
of  10-34.2°C,  while  the  volumetric  heat  capaci¬ 
ties  of  the  various  PCMs  remain  almost  con¬ 
stant.  In  the  theoretical  studies  presented  here 
we  assumed  that  all  phase  transitions  in  the 
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Fig.  5.  Meteorological  data  during  the  experiment  shown  in  Fig.  4. 


PCMs  occurred  in  a  +  1°C  temperature  interval 
around  a  nominal  transition  temperature. 

All  the  parametric  studies  presented  here  were 
performed  using  typical  thermophysical  proper¬ 
ties  for  salt-hydrates,  as  presented  in  Table  1. 

Assuming  that  all  the  collected  daily  solar 
energy  is  stored  as  latent  heat,  the  “melted” 
thickness  of  the  PCM  layer  could  be  estimated 
from  the  simple  equation: 


where:  lpcm  is  the  thickness  of  the  melted  PCM 
layer,  r\  is  the  total  daily  solar  collection  effi¬ 
ciency,  H  is  the  total  daily  solar  radiation  and 
X  is  the  volumetric  latent  heat. 

The  monthly  means  for  the  daily  total  solar 
radiation  in  the  Beer-Sheva  area  ranges  between 
10  MJ/m2  in  January  and  38  MJ/m2  in  June 
(Kudish  and  Ianetz,  1992).  On  the  following 


assumptions:  ( 1 )  in  the  period  for  growing  out- 
of-season  winter  crops  (January-March)  the 
daily  total  solar  radiation  might  reach  a  maxi¬ 
mum  of  25  MJ/m2;  (2)  all  the  solar  energy 
collected  in  the  ISC-PCM  is  stored  as  latent 
heat  of  the  PCM;  (3)  the  volumetric  latent  heat 
capacity  is  255  MJ/m3;  and  (4)  the  solar  collec¬ 
tion  efficiency  ranges  between  65  and  70%,  the 
maximal  thickness  of  that  PCM  layer  that  might 
melt  in  the  winter  growing  season  has  to  be  in 
the  range  of  64-69  mm. 

As  a  case  study  in  this  analysis  we  chose  to 
present  the  results  pertaining  to  the  weather 
conditions  of  January,  which  is  the  coldest 
month.  The  thickness  of  the  PCM  was  chosen 
to  be  67  mm,  as  this  value  was  found  as  first 
estimation  for  the  entire  period  January  to 
March,  for  which  heating  is  required  in  green¬ 
houses.  Meteorological  parameters  were  taken 
from  statistical  data  for  January  in  the  Beer- 
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Fig.  6.  Thermal  energy  storage  per  unit  of  collector  area  vs.  time  of  the  day  for  different  phase  transition 
temperatures. 


TIME,  hour 


Fig.  7.  Location  of  the  interfaces  as  a  function  of  the  time  of  the  day  for  various  phase  transition  temper¬ 
atures.  Tm,  and  Tm>  represent  the  interfaces  between  the  liquid  and  melting  regions  and  between  the  melting 
and  solid  regions,  respectively. 


Sheva  area,  as  presented  in  Fig.  5.  The  overall 
absorption  coefficient  of  the  absorber  metal 
plate  in  January  was  taken  as  0.75,  which  was 
an  average  value  that  was  predicted  from  experi¬ 
mental  measurements  of  the  direct  and  diffuse 
radiation  reflected  from  the  upper  surface  of 
the  ICS  system.  The  convective  heat  transfer 
coefficient  to  the  surroundings  was  estimated  as 
15  W  m-2  K-1. 

The  first  parameter  to  be  studied  was  the 
nominal  melting  temperature.  As  discussed 
above,  the  phase  transition  interfaces  of  the 
PCMs  were  defined  by  the  isotherms  Tm  —  1°C 
and  Tm  +  1°C.  In  terms  of  the  enthalpy 


approach,  an  effective  specific  heat  was  assumed 
to  lie  between  these  two  isotherms. 

The  energy  stored  in  the  ICS  system,  as  a 
function  of  time,  and  the  solar  energy  transferred 
to  the  outer  surface  of  the  ICS  is  shown  in 
Fig.  6.  Figure  7  presents  the  location  of  the 
interfaces  inside  the  PCM.  From  Fig.  6  it  can 
be  seen  that  the  melting  temperature  has  only 
a  minor  influence  on  the  solar  absorption  effi¬ 
ciency.  For  melting  temperatures  in  the  range 
of  15-30°C  the  solar  absorption  efficiency  is 
65-70%.  However,  from  Fig.  7  it  can  be  seen 
that  the  melting  temperatures  have  a  major 
influence  on  the  rates  of  phase  change.  The 
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Fig.  8.  Location  of  the  interfaces  as  a  function  of  the  time  of  the  day  for  various  thicknesses  of  the  PCM 
layer.  The  nominal  phase  transition  temperature  is  20°C.  Tm,  and  Tms  represent  the  interfaces  between  the 
liquid  and  melting  regions  and  between  the  melting  and  solid  regions,  respectively. 


energy  charging  period  may  be  divided  into  two 
stages  (Fig.  7)  as  follows:  The  first  stage  can  be 
defined  as  the  period  in  which  the  two  interfaces 
propagate  in  the  same  direction,  towards  the 
bottom  of  the  container.  This  stage  starts  later 
for  higher  melting  temperature  materials  and 
ends  in  the  afternoon,  approximately  at  15:00. 
The  second  stage  can  be  defined  as  the  period 
in  which  the  solar  radiation  decays  and  the 
temperature  distribution  in  the  PCM  tends  to 
become  uniform  around  its  phase  transition 
temperatures.  In  the  second  stage  the  two  inter¬ 
faces  propagate  in  opposite  directions.  From  an 
examination  of  the  material  melting  at  15°C,  for 
example,  it  can  be  seen  from  Fig.  7  that  by  17:00 
the  solid  front  interface  reaches  the  bottom  of 
the  container  and  by  17:30  the  liquid  front 
interface  returns  to  the  SHTL-PCM  interface. 
Thus,  by  17:30  all  the  PCM  is  within  the  phase 
transition  temperature  interval.  A  numerical 
solution  for  the  same  problem  with  a  melting 
temperature  of  35°C  is  not  presented  here,  since 
for  meteorological  conditions  of  January  and  a 
PCM  thickness  of  67  mm  all  the  energy  is  stored 
as  sensible  heat. 

The  effect  of  the  thickness  of  the  PCM  was 
studied  in  the  range  between  20-100  mm.  We 
chose  a  nominal  melting  temperature  of  20°C 
for  this  parametric  study.  The  locations  of  the 
PCM  interfaces,  as  a  function  of  time  of  the 
day,  for  various  thicknesses  of  the  PCM  are 
shown  in  Fig.  8.  The  maximum  propagation 
distances  that  the  liquid  and  solid  fronts  reached 
during  the  day  for  various  PCM  thicknesses  are 


given  in  Table  2.  The  theoretical  results  showed 
that  for  a  thickness  of  20  mm,  the  liquid  phase 
front  reached  the  bottom  at  about  14:30 
(Table  2),  i.e.  all  the  PCM  had  completely 
melted.  The  maximum  propagation  of  the  liquid 
front  increased  with  decreasing  thickness  of  the 
PCM  layer  (Fig.  8  and  Table  2).  For  a  thickness 
range  of  30-100  mm  this  maximum  value  is 
reached  almost  at  the  same  time  (15:20)  for  all 
cases.  With  regard  to  the  propagation  of  the 
solid  phase  front,  for  thicknesses  less  than 
60  mm  this  front  reached  to  the  bottom  surface 
of  the  PCM  layer,  while  for  thicknesses  above 
70  mm  this  front  did  not  reach  the  bottom  and 
the  maximum  penetration  distance  decreased 
with  increasing  of  PCM  thickness. 

Thus,  under  the  conditions  of  this  study,  the 

Table  2.  Maximal  propagation  distance  of  the  liquid  and 
solid  front  interfaces  during  the  day  for  a  nominal  melting 
temperature  of  20°C  and  average  meteorological  data  for 
January 


PCM  layer 
thickness 

Maximal  phase  front  propagation 

Liquid  interface 

Solid  interface 

Distance 

Time 
of  day 

Distance 

(mm) 

Time 
of  day 

20 

20 

14:30 

20 

10:00 

30 

14.3 

15:20 

30 

11:30 

40 

13.4 

15:20 

40 

13:00 

50 

12.5 

15:20 

50 

15:00 

60 

11.9 

15:20 

60 

17:00 

70 

11.2 

15:20 

45.9 

18:00 

80 

10.7 

15:20 

37.1 

18:00 

too 

10.7 

15:20 

30.0 

18:00 
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thickness  of  the  PCM  layer  had  a  major  influ¬ 
ence  on  the  location  of  the  propagation  distance 
of  the  solid  interface  and  a  minor  influence  on 
the  liquid  interface.  However,  it  should  be  noted 
that  both  the  formation  of  the  interfaces  and 
the  disappearance  of  the  liquid  phase  interfaces 
took  place  at  the  same  time  for  all  four  cases 
examined  in  this  study  (Fig.  8). 

For  the  particular  case  study  of  a  nominal 
transition  temperature  of  20°C  and  meteorologi¬ 
cal  data  for  the  coldest  month  (January),  a 
PCM  layer  thickness  of  no  more  than  30  mm  is 
required  to  store  the  daily  solar  energy. 

5.  CONCLUSIONS 

A  new  type  of  ICS  system  for  a  low-tem¬ 
perature  heat  supply  has  been  presented.  The 
ICS  system  uses  a  salt-hydrate  PCM  as  a  solar 
energy  storage  medium.  The  discharging  process 
is  performed  by  a  finned  coil  heat  exchanger 
located  in  an  organic  heat  transfer  fluid  layer 
(oil)  floating  above  the  PCM.  The  main  advan¬ 
tage  of  this  method  is  that  it  uses  a  water/oil 
heat  exchanger  in  place  of  an  expensive 
water/salt-hydrate  surface  heat  exchanger.  The 
possible  use  of  this  system  is  limited  to  special 
applications,  such  as  space  and  soil  heating  of 
greenhouses  located  in  arid  zone  areas,  during 
winter. 

A  simplified  theoretical  model  to  study  the 
thermal  performance  of  such  a  system  in  the 
charging  process  was  developed,  and  an  exper¬ 
imental  system  was  set  up  to  validate  it. 
Comparison  between  mathematical  results  of 
the  model  with  experimental  data,  for  the  tem¬ 
perature  profiles  in  the  salt-hydrate  PCM  and 
in  the  SHTL,  and  for  the  location  of  the 
solid/liquid  interfaces  in  the  PCM  as  function 
of  time,  showed  reasonable  agreement  for  engi¬ 
neering  design  purposes.  The  comparison  also 
indicated  that  the  model  could  be  improved  by 
taking  into  consideration  the  effect  of  convective 
heat  transfer  in  the  SHTL. 

Our  experimental  study  showed  that  in  the 
winter  in  Beer-Sheva  about  two-thirds  of  the 
total  daily  solar  energy  falling  on  the  collector 
surface  could  be  collected  and  stored  at  28°C 
using  a  PCM  based  on  a  eutectic  mixture  of 
calcium  chloride. 

Theoretical  study  on  the  effects  of  the  trans¬ 
ition  phase  change  temperature  (in  the  range 
15-35°C)  and  the  thickness  of  the  PCM  layer 
(in  the  range  20-100  mm)  showed  that  for  the 
phase  transition  at  20°C — which  is  suitable  for 


heating  greenhouses — the  PCM  thickness  for 
energy  collection  has  to  be  30-65  mm.  However, 
the  optimal  design  of  this  parameter  should  also 
take  into  consideration  the  discharging  process 
and  engineering  aspects  of  the  system;  this  will 
be  undertaken  in  the  future. 
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NOMENCLATURE 

a  effective  absorption  coefficient  of  the  integrated  system 
A  surface  area  of  the  collector  (m2) 
b  transparent  cover  and  air  layer  thickness  (m) 
h  heat  transfer  coefficient  by  convection  (Wm  2  K_1) 
H  total  daily  solar  radiation  (J/m2) 
g  initial  temperature  profile  function 
k  thermal  conductivity  (Wm-1  K-1) 
l  thickness  (m) 

L  the  thickness  of  the  solar  collector  storage  system  (m) 
q  heat  flux  (Wm-2) 

R  thermal  resistance  to  heat  flow  (K.  W  ') 
t  time  (s) 

T  temperature  (K  or°C) 
u  wind  velocity  (ms-1) 
x  coordinate  (m) 
a  thermal  diffusivity  (m2  s"1) 

<5  interface  location  (m) 

X  volumetric  latent  heat  (MJ  m-3) 
r\  total  daily  solar  collection  efficiency 

Indexes 

c  transparent  cover  and  air  layer 
f  stationary  heat  transfer  liquid,  SHTL 
h  convection 
i,  j  layers/sublayer  index 
C  liquid 
m  melting 
r  radiation 

oo  surroundings 
pem  phase  change  material 
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